Late transient acceleration of the universe in string theory on S 1 /Z 2 
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Recently, in Gong et al |T7]] and Wang and Santos 21] it was shown that the effective cosmological 
constant on each of the two orbifold branes can be easily lowered to its current observational value, by 
using the large extra dimensions in the framework of both M-Theory and string theory on S 1 jZ-z. 
In this paper, we study the current acceleration of the universe, using the formulas developed 
in [2J]. We first construct explicitly time-dependent solution to the 10- dimensional bulk of the 
Neveu-Schwarz/Neveu-Schwarz sector, compactified on a 5-dimensional torus. Then, we write down 
the generalized Friedmann equations on each of the two dynamical branes, and fit the models to 
the 182 gold supernova la data and the BAO parameter from SDSS, using both of our MINUIT 
and Monte-Carlo Markov Chain (MCMC) codes. With the best fitting values of the parameters 
involved as initial conditions, we integrate the generalized Friedmann equations numerically and 
find the future evolution of the universe. We find that it depends on the choice of the radion 
potentials V} (I — 1,2) of the branes. In particular, when choosing them to be the Goldberger- 



r{I) 



Wise potentials, V± — \\ tip — Vi J , we find that the current acceleration of the universe driven 
by the effective cosmological constant is only temporary. Due to the effects of the potentials, the 
universe will be finally in its decelerating expansion phase again. We also study the proper distance 
between the two branes, and find that it remains almost constant during the whole future evolution 
of the universe in all the models considered. 

PACS numbers: 03.50.+h, ll.10.Kk, 98.80.Cq, 97.60.-s 
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I. INTRODUCTION 



One of the long-standing problems in particle physics 
and cosmology is the cosmological constant problem: its 
theoretical expectation values from quantum field the- 
ory exceed observational limits by 120 orders of magni- 
tude [1J. Even if such high energies are suppressed by 
supersymmetry, the electroweak corrections are still 56 
orders higher. This problem was further sharpened by re- 
cent observations of supernova (SN) la, which reveal the 
striking discovery that our universe has lately been in its 
accelerated expansion phase I2J]. Cross checks from the 
cosmic microwave background radiation and large scale 
structure all confirm this unexpected result [3| . 

In Einstein's theory of gravity such an expansion can 
be achieved by a tiny positive cosmological constant. In 
fact, such a constant is well consistent with all observa- 
tions carried out so far [J|. Therefore, solving the cosmo- 
logical constant problem now becomes more urgent than 
ever before. As a matter of fact, it is exactly because of 
this that a large number of ambitious projects have been 
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proposed lately to distinguish the cosmological constant 
from dynamical dark energy models [5J. 

Since the problem is intimately related to quantum 
gravity, its solution is expected to come from quantum 
gravity, too. At the present, string/M-Theory is our best 
bet for a consistent quantum theory of gravity, so it is rea- 
sonable to ask what string/M-Theory has to say about 
the cosmological constant. In the string landscape pi], 
it is expected that there are many different vacua with 
different local cosmological constants Q. Using the an- 
thropic principle, one may select the low energy vacuum 
in which we can exist. However, many theorists still hope 
to explain the problem without invoking the existence of 
ourselves in the universe. 

Townsend and Wohlfarth 8] considered a time- 
dependent compactification of pure gravity in higher di- 
mensions with hyperbolic internal space to circumvent 
Gibbons' non-go theorem Q. Their exact solution ex- 
hibits a short period of acceleration. The solution is the 
zero-flux limit of spacelike branes [lfj |. If non-zero flux 
or forms are turned on, a transient acceleration exists 
for both compact internal hyperbolic and flat spaces [ll| . 
Other accelerating solutions by compactifying more com- 
plicated time-dependent internal spaces can be found in 

B 

In the same spirit, the cosmological constant prob- 
lem was also studied in the framework of brane world 
in 5-dimensional spacetimes [13] and 6-dimensional su- 



pergravity [Tj|. However, it turned out that in the 5- 
dimensional case hidden fine-tunings are required [la ]. 
while in the 6-dimensional case it is still not clear whether 
loop corrections can be as small as required [la ]. 

Recently, we [l7[ studied the cosmological constant 
problem and late acceleration of the universe in the 
framework of the Horava-Witten heterotic M-Theory on 
S 1 /Z2 [18]. In particular, using the Arkani-Hamed- 
Dimopoulos-Dvali (ADD) mechanism of large extra di- 
mensions [19( , we showed explicitly that the effective cos- 
mological constant on each of the two orbifold branes can 
be easily lowered to its current observational value. Ap- 
plying it to cosmology, we further found that the domi- 
nation of the effective cosmological constant is only tem- 
porary. Due to the interaction of the bulk and the brane, 
the universe will be in its decelerating expansion phase 
again in the future, whereby all problems connected with 
string cosmology [2(| are resolved. Such studies were also 
generalized to string theory, and found that the ADD 
mechanism can be used in the same way to solve the 
cosmological constant problem [21]. Therefore, the ADD 
mechanism for solving both the cosmological constant 
problem and the hierarchy problem is a built-in mecha- 
nism in the brane world of string/M-Theory. 

In this paper, we apply the theory developed in [211 ] 
in the framework of string theory on S 1 jZ^ to cosmol- 
ogy, and study the current acceleration of the universe. 
In particular, in Sec. II we first give a brief review 
of the theory, and then write down the field equations 
both in the bulk and on the two orbifold branes. In 
Sec. Ill, we present a particular time-dependent solution 
to these equations, and study its local and global prop- 
erties. Then, we write down explicitly the generalized 
Friedmann equations on each of the two branes for any 
radion potentials V± {I ■= 1,2) of the branes. Depend- 
ing on the choice of V^ s, the future evolution of the 



universe is different. We study two different cases. We 
fit these models to the 182 gold supernova la data [22j 
and BAO parameter from SDSS J23(, using both of our 
MINUIT [H and Monte-Carlo Markov Chain (MCMC) 
[251 ] codes. With these best fitting values as the initial 
condition, we integrate numerically the field equations on 
the branes to find the future evolution of the universe. In 
the latter case, we show that the current acceleration of 
the universe driven by the effective cosmological constant 
is only temporary. Due to the effects of the potentials, 
the universe will be in its decelerating expansion phase 
again. We also study the proper distance between the 
two orbifold branes, and find that it remains almost con- 
stant during the whole future evolution of the universe 
in all these models. In Sec. IV, we summarize our main 
results and present some concluding remarks. 

Before turning to the next section, we would like to 
note that Sahni and Shtanov [20] found that transient 
acceleration of the universe happened in the DGP brane 
model [27J], too. 



II. BRANY COSMOLOGY OF STRING 
THEORY ON S 1 /Z 2 

To begin with, in this section we give a brief review on 
the cosmological models of orbifold branes developed in 
[211 ] . In this paper we shall restrict ourselves directly to 
the case D = d = 5. For the case with arbitrary D and 
d, we refer readers to [211 . 



A. The Model 

For the toroidal compactification of the Neveu- 
Schwarz/Neveu-Schwarz (NS-NS) sector in (5+5) dimen- 
sions, Mio = M5 x 75, where 7r is a 5-dimensional torus, 
the action takes the form 28, 23, 
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V A $) (V B $) - ^i? 2 | , (2.1) 



where V a denotes the covariant derivative with respect 
to g AB with A, B = 0, 1, „., 9, and $ is the dilaton field. 
The NS three-form field Habc is defined as 

Habc — 3d[ A B BC ], (2.2) 

where the square brackets imply total antisymmetriza- 
tion over all indices. The 10-dimcnsional spacetimes to 
be considered are described by the metric, 

^10 = 9Asdx dx 

= g ab (x c )dx a dx b + h tj (x c )dz z dz j , (2.3) 

where g ab is the metric on M 5 , parametrized by the co- 
ordinates x a with a, b, c = 0, 1, ...,4, and hij is the metric 
on the compact space T§ with periodic coordinates z l , 
where i, j — 5,6, ..., 9. 

By assuming that all the matter fields are functions of 
x a only, it can be shown that the effective 5— dimensional 
action is given by, 

S 5 = -^Jd 5 x^/\g7\e-^{R 5 [g} + (^a4>){^ a 4>) 

+ \ (V a h ij ) (v°fc«) - ^H ab cH a 

VaBij) (\7 a B kl 
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with the 5— dimensional internal volume given by 

V(x a )= fd 5 z^\h\ = \h\ 1 ^V Q . 



(2.4) 

(2.5) 
(2.6) 

(2.7) 



Note that in writing the action (|2.4|) we had assumed 
that the flux is block diagonal, 



(b cd ) 
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(2.8) 



The action (|2.4[) is usually referred to as written in the 
string frame. To go to the Einstein frame, we make the 
following conformal transformations, 
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Then, the action (12.41) takes the form 
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(2.10) 



where V a denotes the covariant derivative with respect 
to g ao . It should be noted that, since the definition of 
the three-form Habc given by (|2.2j) is independent of 
the metric, it is conformally invariant. In particular, we 
have H a bc = H aoc and B a b = B ao . However, we do have 



H" 



ad be cf tj c\ — 6 ijabc 

g g g H def = il H , 



tj ijabc r) — 6 tj ijabc 



(2.11) 



Considering the addition of a potential term [29J, in 
the string frame we have 



STo = - I d 10 xy/\§w\Vf . (2.12) 



Then, after the dimensional reduction we find 

S 5 , m = -Vo J d 5 xy%\ \h\ 1/2 V&, (2.13) 

where 
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h 
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we find that 
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where the effective 5— dimensional potential (|2.14j) now 
becomes 



V 5 = V ( ° 5) exp 



V6 



i> 



(2.19) 



where V^ = 2k 2 5 V V 1 s . 

To study orbifold branes, we consider the brane ac- 
tions, 
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(2.20) 



where J = 1,2, V4 (<fi,ip) denotes the potential of the 
scalar fields <f> and ip, and Cm' s are the intrinsic co- 
ordinates of the I-th brane with /j,, v = 0,1,2,3, and 
(-1 = — e 2 = 1. x denotes collectively the matter fields, 
and g s is a constant, which is related to the four- 
dimensional Newtonian constant via the relation given 
by Eq. (|2.39|) below. The variation of the total action, 
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(2.14) 



Changed to the Einstein frame, the action (|2.13p finally 
takes the form, 
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with respect to the metric g a b yields the field equations, 
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where S(x) denotes the Dirac delta function normalized 
in the sense of [30] , and the two branes are localized on 
the surfaces, 



$/ (x a ) = 0. 



(2.23) 
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The energy-momentum tensors T^ b and T^ v are given 
by 
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and g/ti/ is the reduced metric on the I-th brane, defined 
as 
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Variation of the total action Eq. (|2.21[) with respect to 
</>, ip and B, respectively, yields the following equations 
of the matter fields, 
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<z6 i,v t>~ SB ab 

To write down the field equations on the branes, one 

(5) 

can first express the delta function part of G ab in terms 
of the discontinuities of the first derivatives of the metric 
coefficients, and then equal the delta function parts of the 
two sides of Eq. (|2.22p . as shown systematically in [3l0 | . 
The other way is to use the Gauss-Codacci equations to 
write the (4)-dimensional Einstein tensor as [32j ]. 



G (4) = 0(6) , E {5) , jr(4) 
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where rt a denotes the normal vector to the brane, G'- 5 -' = 

, (5) (5) 

9 G ab , and C^ 6 ;l d the Weyl tensor. The extrinsic curva- 
ture K^ v is defined as 



e? u ^V n 6 . 



(2.34) 



A crucial step of this approach is the Lanczos equations 
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where 
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= 5 ( 7 )^ 






(2.36) 



Assuming that the branes have Zi symmetry, we can 
express the intrinsic curvatures K^J in terms of the ef- 
fective energy-momentum tensor 7^„ through the Lanc- 

zos equations ()2.35j) . Then, we find that G}J given by 
Eq. (|2.32[) can be cast in the form, 

where 

i 1 



(2.37) 
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For a perfect fluid, 



< /ii/ 



(p + p)u fl u„-pg^, 



(2.39) 



(2.40) 



where m m is the four- velocity of the fluid, we find that 



{p + p)u^u v 



P+2 p } ''''" 



(2.41) 



Note that in writing Eqs. (|2.37ll - (|2.4ip . without causing 
any confusion, we had dropped the super indices (7). 

In the rest of this paper, we shall turn off the flux, i.e., 
Bcd — 0, which is consistent with the field equations, 



provided that &),' = and $ 



i(') 



0. 



B. The General Metric of the Five-Dimensional 
Spacetimes 

Since we shall apply such spacetimes to cosmology, let 
us first consider the embedding of a 3-dimensional spatial 
space that is homogeneous, isotropic, and independent of 
time. It is not difficult to show that such a space must 
have a constant curvature and its metric takes the form 

0, 



,^9 dr 2 

dY? k = — 

h 1 - kr 2 



r 2 (d9 2 + sin 2 



') 



(2.42) 



where the constant k represents the curvature of the 3- 
space, and can be positive, negative or zero. Without 
loss of generality, we shall choose coordinates such that 
k = 0, ±1. Then, one can see that the most general 
metric for the five-dimensional spacetime must take the 
form, 

dsl = g ab dx a dx b = g M Ndx M dx N ~ eM*™)^ 2 , (2.43) 



where M, N = 0, 1. Clearly, the metric (|2.43j) is invariant 
under the coordinate transformations, 



jN _ rN (M 



= f N (x M )- 



(2.44) 



Using these two degrees of freedom, without loss of gen- 
erality, we can always set 

3oo =5ii: .9oi = 0, (2.45) 

so that the five-dimensional metric finally takes the form, 

ds 2 = e 2 < 1 ^ (dt 2 - dy 2 ) - e 2 ^^d£ 2 . (2.46) 



It should be noted that metric (12.46|) is still subjected to 
the gauge freedom, 

* = f(t' + v') + g(t' - v% v = /(*' + v') - s(f - v% 

(2.47) 
where f(t' + y') and g(t' — y') are arbitrary functions of 
their indicated arguments. 

It is also interesting to note that in [31j a different 
gauge was used. Instead of setting g 00 — g n it was chosen 
that the two branes are comoving with the coordinates, 
so that they are located on two fixed hypersurfaces y = 
0, y c . For details, see |31j . 



C. The Field Equations Outside the Two Orbifold 
Branes 



The non- vanishing components of the Ricci tensor out- 
side of the two branes are given by 



(5) 

tt 

(5) 



a tyy + 3<T, v 0J !y - [<r iU + 3w jtt + 3w t (w )t - a t )] , 



ty ~ ^3 [u tty + W tt U,y ~ {cr,t^,y + ffjW.t)] , 

R yy = °> + 3cr ,*^t _ Kw + ^ u ,vv + 3w ,a ( w ,y _ a ,v)\ > 
Ri mli = -e~ 2<7 g mn {u,tt + 3w t 2 - (u! tVV + 3u t y 2 ) 

+2ke 2{a - uj) \ , (2.48) 

where now m, n = r, 9, (p, a.t = da / dt and so on. Then, 
it can be shown that outside of the two branes the field 
equations have four independent components, which can 
be cast into the form, 

^,tt + W,t (W t - 2a }t ) + U.yy + LU^y (tjj y y ~ 2(J,y) 

= -\ [(</>/ + </>,y 2 ) + (4>/ + ty 2 )] , (2-49) 



2a, u + w,« ~ 3w, 4 2 - (2a tVV + Lu, yy - 3w y 2 ) - 4fce 2(<T - w) 
= -\[^,t 2 -^y 2 ) + {^,t 2 -^,y 2 )}, (2.50) 

w itt + 3u>/ - (u m + 3u t y 2 ) + 2ke 2{a - u)) 
3 5 ' 



(2.51) 



(2.52) 



where Vs is given by Eq. (|2.19j) . On the other hand, the 
Klein-Gordon equations (|2.27[) and (|2.28p outside the two 
branes take the form, 



4>,tt + 3(/) t tU,t - (4>,yy + 34>,yU, y ) 
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Vse", 



Ip.tt + 3lp,tU,t - {lp,yy + 3lp t yUJ t y) 



-\I^V,e 2 °, 



(2.53) 



(2.54) 



D. The Field Equations on the Two Orbifold 
Branes 



Eqs. (|2.48[) - (|2.54|) are the field equations that are valid 
in between the two orbifold branes, 3/2(^2) < V < 2/i(£i)) 
where y = yi(tj) denote the locations of the two branes. 
The proper distance between the two branes is given by, 



V 



y/-g yy dy. 



(2.55) 



On each of the two branes, the metric reduces to 

ds 2 5 \ M m = fl&tefodefo = drf - a 2 (77) cE 2 , (2.56) 

where £?p. = {77, r, 6, ip}, and 77 denotes the proper time 
of the I-th brane, defined by 



dri 



dtr 



a(ri) = expjw [£7(77), 2/7(77)]} 



(2.57) 



with yj = dyj/dri, etc. For the sake of simplicity and 
without causing any confusion, from now on we shall 
drop all the indices "I", unless some specific attention 
is needed. Then, the normal vector n a and the tangen- 
tial vectors e? •, are given, respectively, by 



n a 
, a 

o a 

-(e) 



ee^i-ySl+iSy), 
-e(y6? + i6 a y ), 
i$t+y s yi e a {r) =5 a r , 



'(¥>) 
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where e = ±1. When e = +1, the normal vector n a 
points toward the increasing direction of y, and when 
e = — 1, it points toward the decreasing direction of y. 
Then, the four-dimensional field equations on each of the 
two branes take the form, 



H z 



^-^(p+n^ + ^ + yp + ti* 
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where _ff = d/a, pa = A4/(87rG4), and 
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with </>.„ = n°V a 0. If the typical size of the extra dimen- 
sions is R, then it can be shown that 



A 4 {R\ w (M 1Q \ W M 4 



(2.62) 



where M p i and Z p / denote the Planck mass and length, 
respectively. If M\q is in the order of TeV [35| , we find 
that, in order to have p\ be in the order of its current 
observations value pa — 10~ 47 GeV 4 , the typical size of 
the extra dimensions should be R ~ 10~ 22 to, which is 
well below the current experimental limit of the extra 
dimensions [36j. 



III. A PARTICULAR MODEL 

In this section, we consider a specific solution of the 
five-dimensional bulk and the corresponding Friedmann 
equations on the orbifold branes. 

A. Exact Solutions in the Bulk 

It can be shown that the following solution satisfies the 
field equations in the bulk, 

««) = iln(t) + iln(£), 



A 



where l>o = \/3/2 4>q + \Jhj2 ipo , from which we find 




t = 

FIG. 1: The Penrose diagram for the metric given by Eq, (|3.3[) 
in the text, where the spacetime is singular at t = 0. The 
curves OP A and OQA describes the history of the two orb- 
ifold branes located on the surfaces y — yi(ri) with I — 1,2. 
The bulk is the region between these two lines. 



10 , , x 
w (t) = _ ln(t), 



for fc = — 1, where 



-£K j 



ln(t)+Vo, 



>n:'M V^ 



(5). 



(3.1) 



(3.2) 



with V'o being an arbitrary constant. Then, the corre- 
sponding 5-dimensional metric takes the form, 

ds 2 = (V\ t 2 ' 9 (dt 2 - dy 2 ) - t^dX 2 ^. (3.3) 

Clearly, the spacetime is singular at t = where all the 
four spatial dimensions collapse into a point singularity, 
a big bang like. This can be seen more clearly from the 
expression, 



ta^ = f<M' Q = |t- 2 °/ 9 . 



(3.4) 



The corresponding Penrose diagram is given by Fig. 1. 

Lifting the solution to the 10-dimensional superstring 
spacetime, we find that in the string frame the metric 
takes the form, 



°^io = gAsdx dx 



= e VF *o | QA t -l/3 ^2 _ dy 2) _ i 5/3 d ^ i 

-&ft* a t~ 1/ *8 i jdz i dz i . (3.5) 



The corresponding dilaton field is given by 

5 
fc = --ln(t) + fc , 



(3.6) 



* *A 50 



0O+-5/3 



-j?>o t 



(3.7) 



Clearly, it is also singular at t = 0, but with a weaker 
strength in comparing to that of the five-dimensional 
spacetime given by Eq. (|2.19| ). A critical difference is 
that in the string frame the proper distance along the 
y-direction becomes decreasing as t increases, in contrast 
to that in the Einstein frame, as can be seen clearly from 
Eqs.O and (|3"3]1 . 



B. Generalized Friedmann Equations on The 
Branes 



On the other hand, from Eq. (|2.61[) we find that 



E& = - 



1 



Qf ] - 



42a 2 
20 



g& = 



31 



126a 2 ' 
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81a 9 / 5 y 378a 2 ' 
where now a(r) = t 10 ' 9 (r), and y is given by 



U 



9/10 



H z - 



6 

7^ 



1/2 



(3.8) 



(3.9) 



with e y = ±1. Inserting Eqs. (|3.8[ ) and (|3.9p into 
Eqs. (|2~l)9"j) and (|2~TJ0"]) . we find that 

H 2 



8wG , , 200 

0> + W)+PA)+7g^2 



3 

+ 3p7^ +T ^) • 

— r— (opA + •J'>~(0,V>) ■ 



(3.10) 



2p - 5 P ) 



2ttG 
3/Oa 



(p + T(<M>)) +P(2p + 3p) 



+ (p + 3p-r (0>v , ) )r( 0:V , ) ] 



(3.11) 



It is remarkable to note that these two equations do not 
depend on both e defined in Eq. (|2.58p and e y defined in 
Eq.flnU). Combining Eqs. (|3~TU|) and (|3~TTj) . we obtain 



{f> + T {4>/l p)) +3H(p + p) 



H 

20A 

(p + 



[4 (p + PA 



~(<t>,i>)l 



PA 



where 



A = ! + ^ (P + T (4>,i>)) 



(3.12) 
(3.13) 



Eq. (|3.12p shows clearly the interaction among the matter 
fields confined on the branes and the bulk. This can also 
be seen from Eq. 



Current Acceleration of the Universe 



To study current acceleration of the universe, we first 



set 

p = 0, 

and then introduce the quantities, 
a - En a - T (^) 

Per Per 



(3.14) 



A = ^, 

Per 



n k = 



200 



n 



(") 



189# 2 a 2 



(3.15) 



where p cr = 3Hq/8ttG. It should be noted the slight 
difference between Sl k defined here and the one normally 
used, Sl k = -k/(H$a 2 ). Then, Eqs. (|3~T0l) . (j3~T2j) and 
(13.91) can be written as 



e 2 = n A + n t + a k 



Sit 

4ft A ; 



sk 



+ 



- - 

A\5 

Sit 

20Sl A 



(0 A + 16fi t - 15fi T ) 
(31Q t - 30Sl T ) 

9/20 



(3.16) 



(3.17) 




x\ n A + n t + 



n. 



40 A ' 



where 25 = H/H Q , y* = dy/d(H T), and 
with the constraint, 



1 = Sl 



(o) 



si A + si\ 



(0) , ^t 



(oy 



40a 



(3.18) 
(3.19) 

(3.20) 



)(0) 



where fl N s denote their current values. On the other- 
hand, in terms of Si's, we find 

— = A (3H A - 2Sl t + 5Sl T ) 
a 10 



~40ft A 



(70 t - 10Sl T ) . 



(3.21) 



To study Eqs. ([3~T6| - ([3~T8l) and (j3~2l"T) further, we need 
to specify Sl T . In the following, we shall consider two 
different cases. 



1. V^=V&exp{f(^tf+v1>)} 

If we choose the potential V^ (0, ip) on each of the two 
branes as [cf. Eq. (f2~T9|) ]. 



"i"^Wl ^ 




(3.22) 




FIG. 2: The marginalized contour of Q m — Qa for the potential 
given by Eq. (|3.22[l with n = 1. 



where V, A \ an d n are arbitrary constants, we find that 



'(4) 



f7 r=e/ ^ 



n/2 



n 



(0) 



*%. 



(3.23) 



Then, our fitting parameters in this case can be chosen 
as 



{^«4 0) } 



(3.24) 



for any given n. 

Fit ting the above model to the 182 gold supernova la 
data [22| and the B AO parameter from SDSS [23| , by us- 
ing our numerical code [2J|, based on the publicly avail- 
able MINUIT program of CERN, we find that, for n = 1, 
the best fitting is Sl m = 0.24±ggg, Sl A = 0.76 ±g|?, and 
fi/c = 0.00 ± o 1 ;^ with x 2 = 172.4. Figs. 2-4 show the 
marginalized contours of the Si's, from which we can see 
that the effect of the interaction between the bulk and 
the brane is negligible, and the later evolution of the uni- 
verse follows more or less the same pattern as that of the 
ACDM model in the Einstein theory of gravity. 

For n — 3.5, we find that the best fitting is Sl m = 
0.27 ± g;g|, Sl A = 0.58 ± ° H, and Sl k = 0.00 ± g : gg with 
X 2 = 164.2. Figs. 5-7 show the marginalized contours of 
the Si's. 

The above shows clearly that the case with n = 3.5 is 
observationally more favorable than that of n = 1. We 
have also fitted the data with various values of n, and 
found that the best fitting value of n is about n = 3.5. 

With the above best fitting values of the Si's and n as 
initial conditions, the future evolution of the universe is 
shown in Figs. 8 and 9, from which we can see that all 
of them, except for SI a, decreases rapidly, and Sl A soon 
dominates the evolution of the universe, whereby a de 
Sitter universe is resulted. 

From the metrics of Eqs. (|3.3p and (|3.5p . on the other 
hand, one may naively conclude that the radion in the 





FIG. 3: The marginalized contour of f2 m — fi*. for the potential 
given by Eq. ((3~2"2"j) with n = 1. 



FIG. 6: The marginalized contour of f2 m — fi*. for the potential 
given by Eq. ((3~2"2"T) with n = 3.5. 





FIG. 4: The marginalized contour of Qk — ^a for the potential 
given by Eq. (pQ2J) with n = 1. 



FIG. 7: The marginalized contour of fife — Qa for the potential 
given by Eq. (|3.22|l with n = 3.5. 




0.26 0.28 



0.3 0.32 0.34 



x 

a 




1 1.2 1.4 1.6 U 

a/a„ 



2.2 2.4 



FIG. 5: The marginalized contour of fi m — J7a for the potential 
given by Eo. (|3.22p with n = 3.5. 



FIG. 8: The evolution of the matter components, £Vs, for the 
potential given by Eq. (|3.22[) with n = 3.5. 
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0.22 0.26 0.3 



0.5 0.6 0.7 



-0.2 0.2 0.4 
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FIG. 9: The evolution of the acceleration a** /a = 
(d 2 a/d(HoT) )/a for the potential given by Eq. (|3.22|l with 
n = 3.5. 



FIG. 11: The marginalized probabilities and contours for the 
potential given by Eq. (|3.22[) with n = 3.5. 
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FIG. 10: The locations of the two branes yi(a), and the proper 
distance T> between the two branes for the potential given by 
Eq. (|3.22|l with n = 3.5. The initial conditions are chosen so 
that yi(ao) = 3 and 3/2(00) = 1. The choice of e y = +1 
(e y = — 1) corresponds to the case where the branes move 
towards the increasing (decreasing) direction of y. 



present case is not stable, as the proper distance given by 
Eq. (|2.55[) seems either to increases to infinity (in the Ein- 
stein frame, given by Eq. (|3.3p ) or to decreases to zero (in 
the string frame, given by Eq. (|3.5p ). as t — » 00. A closer 
investigation shows that the problem is not as simple as 
it looks like. In particular, since yj = yi(ri), Eq. (|2.55|) 
makes sense only when the relation n = Ti(t^) is known. 
In the present case, we transform such a dependence to 
the expansion factor a, and plot it out in Fig. [TU1 together 
with yi(a), from which we can see clearly that the dis- 
tance between the two branes remains almost constant. 
This indicates that the radion might be stable. Certainly, 
before a definitive conclusion is reached, more detailed 
investigations are needed. 



We also fit the above model with n = 3.5 _by using our 
Monte-Carlo Markov Chain (MCMC) code J2{|, based on 
the publicly available package COSMOMC [371] , and find 
that the best fitting is ft m = 0.27±22t ft 
and ft,. = -0.0026 ± 



-0.03' 



n ci -1-0.09 
U - Di=n 0.10' 



alio with X 2 = 164.10, where 



ft* 



ft?- 



(3.25) 



The corresponding marginalized probabilities and con- 
tours are given in Fig. [TTJ Clearly, these best fitting 
values are consistent with those obtained above by using 
our MINUET code (2I. 



2. vy^wtf-vsy 

To stabilize the radion, Goldberger and Wise proposed 



r(I) 



to choose the potential V4 as [38] , 



vfW) = Ai /) (V a -V) 2 



(3.26) 



where A4 and vi 2 (I = 1,2) are constants. Then, we 
find that 



fil 7 ) = ft!°< 7 ) 




ln(a) -vj 2 



(3.27) 



where fti- = ei\ 4 / p cr - Note that in writing the above 
expressions, without loss of any generality, we had set 
ip(s = 0. Then, the fitting parameters now can be taken 

as. 



{ft A ,ft(°>,fti°U} 



(3.28) 



Fittingthe above model to the 182 gold supernova la 
data [22| and the BAO parameter from SDSS [23|, we 
first study the dependence of x 2 on v i- Table U shows 
such a dependence and the best fitting values of ftj's for 
each given vi. 
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VI 


x 2 


i 'm 


n k 


«A 


Qa + T 


0.1 


171.28 


0.25=b ; 04 


-0.0009±g;i 


o.72±8:81 


0.73 


0.3 


168.10 


o.29±8:Sl 


-o.ooo6±8:S 


l-06±g;l? 


0.47 


0.5 


157.50 


0.29±g;g| 


-0.008±[ji| 


1.28± ; 28 


0.70 


1.0 


156.69 


o.29±8:o4 


-o.oo2±8:ll 


i.64±8:2 


0.64 


3.0 


156.38 


o.28±8:o4 


-0.008±g;I| 


1.93±J$ 


0.57 


10.0 


166.35 


o.28±8:Sl 


-o.oo2±8:ll 


1.97±§# 


0.56 



TABLE I: The best fitting values of Qi for a given vi of the 
potential given by Eq. (|3.26[) . 




FIG. 12: The marginalized probabilities and contours for the 
potential given by Eg. (|3. 26(1 with vi = 0.5. 



From the table we can see that \ 2 decreases until 
vi ~ 3.0 and then starts to increase, as vi is continu- 
ously increasing. However, fl\ and its uncertainty also 
increase as vj is increasing, while Q m and Qk remain al- 
most the same. Since f2 T acts as a varying cosmological 
constant, Table [I] shows that the total effective cosmo- 
logical constant £1^ = Qa + ^t is between 0.47 and 
0.73. 

Fig. [T2] shows the marginalized probabilities and con- 
tours for the potential given by Eq. (|3.26p with vi = 0.5, 
and Fig. [TBI shows the future evolution of the correspond- 
ing acceleration of the universe. From there we can see 
that the acceleration increases to a maximal value and 
then starts to decrease. As the time is continuously in- 
creasing, it will pass the zero point and then becomes 
negative. Thus, in the present model, the domination of 
the cosmological constant is only temporary. Due to the 
presence of the potential term, represented by f2 r , the 
universe will be in its decelerating expansion phase again 
in the future, whereby all problems connected with a far 
future de Sitter universe are resolved [20]. The effects of 
f2 T can be seen clearly from Fig. 1141 from which we can 
see that both O m and flk decrease rapidly, and soon Q T 
dominates the evolution of the universe. 

These are the common features for any given value of 
vj. Figs. 11511161 and 1 171 show, respectively, the marginal- 




FIG. 13: The acceleration a** /a for the potential given by 
Eq.fl32SJ) with V! = 0.5. 
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FIG. 14: The future evolution of Qi for the potential given 
by Eq.j3]26) with Vi = 0.5. 



ized probabilities and contours, the future evolution of 
a** J a and of Oj for vj = 0.1. 

In addition, we also find that the proper distance be- 
tween the two orbifold branes defined by Eq. (|2.55p is not 
sensitive to the choice of vi, and remains almost constant 
during the future evolution of the universe, as shown in 
Fig. [T8l This also indicates that the radion might be 
stable in the present case, too. 



IV. CONCLUSIONS AND REMARKS 

Recently, we studied the cosmological constant prob- 
lem in the framework of both M-theory 17J and string 
theory [2l| on S 1 /Z2, and showed that, among other 
things, the effective cosmological constant on the branes 
can be easily lowered to its current observational value 
using the ADD large extra dimension mechanism [19j |. 
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FIG. 15: The marginalized probabilities and contours for the 
potential given by Eq. (|3.26p with vi = 0.1. 
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FIG. 16: The acceleration a"* /a for the potential given by 
Eq. p^o) with vi = 0.1. 
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FIG. 17: The future evolution of Qi for the potential given 
by Eq.([3j26) with vi = 0.1. 
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FIG. 18: The locations of the two branes, yi(a), and the 
proper distance, T>, between the two branes for the potential 
given by Eq. (|3.26[) with Vi = 0.5. The initial conditions are 
chosen so that yi(ao) = 3 and 3/2(00) = 1. The choice of 
e y — +1 (e y — —1) corresponds to the case where the branes 
move towards the increasing (decreasing) direction of y. 



Thus, brany cosmology of string/M- Theory seems to 
have a built-in mechanism for solving both the cosmo- 
logical constant and the hierarchy problems. 

In this paper, we have studied a particular cosmolog- 
ical model in the framework of string theory on S 1 /Z 2 , 
developed in [2l| . We have first solved the field equations 
in the bulk and then studied its local and global proper- 
ties. In particular, we have found that the 10-dimensional 
bulk has a big-bang-like singularity at t = 0. 

After obtained explicitly the generalized Fried- 
mann equations on each of the two branes for 



any radion potentials V 4 

branes, we have studied 

( 7 ) T/0 Qv ^ / n ( 5 



two 



= 1 



different 



2) 
cases 



of the 
where 



V. 



V°4) ex P 



2 \V6 



*)} 



and V A 



(i) 



W hP 2 



VI 



respectively. For each of these poten- 



tials, we have fit the corresponding models to the 182 
gold supernova la data [22| and the BAO parameter from 
SDSS [23|, and obtained the best fitting values of the pa- 
rameters involved. To doubly check our numerical codes, 
we have used both of our MINUIT [H and Monte-Carlo 
Markov Chain (MCMC) [25| codes, and gotten the same 
results within the allowed errors. With these best fit- 
ting values as the initial condition, we have integrated 
numerically the field equations on the branes, and found 
the future evolution of the universe. In particular, for 



V A 



(J) 



\ ( i } (v> 2 



■vr 



we have found that the current 



acceleration of the universe driven by the effective cosmo- 
logical constant is only temporary. Due to the effects of 
the potentials, the universe will be in its decelerating ex- 
pansion phase again in the future. We have also studied 
the proper distance between the two orbifold branes, and 
found that it remains almost constant during the whole 
future evolution of the universe in all these models. 
In the framework of orbifold branes, an important 



question is the radion stability. The considerations of 
the proper distance between the two orbifold branes in- 
dicate that in these cases the radion might be stable, 
although further studies are highly demanded. Recently, 
two of the authors (NOS & AW) studied the problem 
in a static background with a four-dimensional Poincare 
symmetry, 



where 



m = 

V>(y) = 



e 2aM (v^dx^dx" - dy 2 ) , 



(4.1) 



-ln|y + t/o|, 



:ln|y + y | + <fo, 



— ln\y + y \ + ip , 



(4.2) 



and j/Oj 4>o & n d V'o are the integration constants with, 



V"o 



In 



9F° 

(5). 



7E" 



(4.3) 
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Following [39(, we are currently investigating the radion 
stability, and wish to report our findings some where else 
soon. 
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